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ABSTRACT 

IZwl8, a blue compact dwarf (BCD) galaxy with the 2nd lowest metallicity measured in a star- 
forming object, has been observed with all three instruments on board the Spitzer Space Telescope. 
We present the deepest 5—36 /im mid-infrared (mid-IR) spectrum of this galaxy as yet obtained, as well 
as 3.6 /xm to 70 /im imaging results. As with SBS 0335-052E, another BCD with similar metallicity, 
IZwl8 shows no detectable emission from polycyclic aromatic hydrocarbons (PAHs). However, the 
continuum emission, from 15 to 70 /im, of I Zw 18 has a much steeper slope, more characteristic of 
a typical starburst galaxy of solar abundance. The neon abundance as measured from the infrared 
fine-structure lines is ~ 1/23 Z Q , and the sulfur abundance is ~ 1/35 Z Q , generally consistent with the 
nebular oxygen abundance of 1/30 Zq derived from optical lines. This suggests that the extinction 
to the infrared emitting regions of this galaxy is low, also in agreement with the optical Balmer line 
ratios. 

Subject headings: dust, extinction — galaxies: individual (IZwl8) — galaxies: dwarf — galaxies: 
starburst 



1. INTRODUCTION 

Low-metallicity galaxies may have been the first sites of 
star formation in th e early universe ([White fc ReedH97§ ; 
iDekel fc Silld [l986). Consequently understanding their 
properties may provide valuable insights into the forma- 
tion of the first generation of normal stars and the subse- 
quent enrichment of the interstellar medium. However, 
finding truly primordial galaxies in the distant universe is 
beyond the reach of current technology. At high redshift, 
observational limitations introduce biases toward the de- 
tection of high-mass high-luminosity systems in which 
the short time-scales of massive star formation lead to 
the identificati on of systems that a re already chemically 
enriched (i.e., Maiolino et al.ll2~003ft . An alternative ap- 
proach is to identify and study unevolved galaxies in the 
local universe. Such a sample are the blue compact dwarf 
galaxies (Isee the review bv Kunth fc O stlin 200(3). 

Since its discov ery bv IZwickvl (1 19661) and the 
seminal paper of Sca rle fc Sargent] (11972ft . IZwl8 
has been studied extensively at many wavelengths. 
With an oxygen abundance determined from the 
optical lines in Hii regions of 12+log(Q /H)=7.17 



/H) 

(jSkillman fc Kennicutt|[l993t llzotov fc Thuanl 11999ft . or 

~1/30Z Q 6 , it has remained the lowest metallicity BCD 
for over two decades until the recent study of the 
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western component of SBS 0335-052 (jlzotov et al J 120051 : 
Papadcros ct al. 2006) , which has slightly lower metallic- 
ity (12+log(Q/H)=7. 12). Distance estimates to IZwl8 
range from ~10 Mpc (Hun teTfc Thronsonlll995ft . ~12.6 
Mpc (|Ostliijl2000ft . and up to ~15 Mpc (jlzotov fc Thuanl 
12004ft . Here we adopt a distance of 12.6 Mpc (1"~61 pc). 
IZwl8 consists of two bright knots of star-formation, 
a northwest (NW) component and a southeast (SE) 
one, together they form the "main body" of the sys- 
tem. Both the NW and SE components contain nu- 
merous young star clusters, with ages ranging from 3 
to lOMyr, but the age of the underlying stellar popu- 
lation in I Zw 18 is still a matter of debate , with max- 
imum ages ranging from 500 Myr to 5 Gyr (|Aloisi et al.l 
[19991 lOstlinll2000t iRecchi et al J 12001 [Hunt et al.l l20ol 
llzotov fc Thuanl 120041 : IMomanv et~atl 12005ft . Approxi- 
mately 22" northwest of the "main body" , there is a 
blue irregular low surface brightness star-forming region, 
called the "C com ponent" , which is embedded in a com- 
mon Hi envelope (jvan Zee et alJ fl998t llzotov fc Thuanl 
l200l . 

It is well known that some BCDs are forming stars 
at a rate that can only be maintained for ~1/10 
of Hubble time given their store of available hydro- 
gen. SBS 0335-052E, which has a similar metallicity 
(12+log[0/H|=7.31), is forming star s at a high rate in 



Crete, 



2 University of 
Heraklion, Greece 

3 IESL/Foundation for Research and Technology 



Hellas, GR- 
Observatoire do 



71110, Heraklion, Greece, and Chercheur Associe 
Paris, F-75014, Paris, France 

4 INAF-Istituto di Radioastronomia, Largo E. Fermi 5, 50125 
Florence, Italy 

5 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 
RA Leiden, The Netherlands 

6 llzotov 



Thuanl 0999) give 1/50 Z 



©• 



Here, 
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the new oxygen solar abundance of 12+log[0/H]=8.69 
IjAllend c Prieto ct al. 2001), which results in a metallicity of 



1/30 Z 
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forming stars at a slo wer rate in compl exes which are 
diffuse and extended (|Hunt et alj|200 5a). These differ- 
ences have been attributed to different star-formation 
modes. SBS 0335-052E is in an "active" mode with rela- 
tively high star-formation rate (SFR) in compact dense 
regions, while in I Zw 18 star-formation occurs in a "pas- 
sive" mo de in more extended reg ions with a relatively 
low SFR (|Hirashita fc Huntj|200l . 

In this paper, we present a detailed analysis of the in- 
frared properties of the "main body" of I Zw 18 based on 
the current deepest mid-IR spectra for this galaxy, ob- 
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tained using the infrared spectrograph (TRS .lHouck et al~ 
2004al) on the Spitzer Space Telescope ([Werner et al. 
2004) . We also discuss the results of pho tometric ob- 
servations with I RAC (|Fazio et all I2004D and MIPS 
(jRieke et al.l l2004f ) observations. We describe the ob- 
servation and data analysis in §2. The Spitzer images 
and spectra are discussed in §3, along with the observed 
morphologies and abundance estimates. A discussion of 
the SFR of I Zw 18 is given in §4 and we summarize our 
conclusions in §5. 

2. OBSERVATIONS 
2.1. Spitzer /IRS Spectroscopy 

I Zw 18 was observed as part of the IRS 7 Guaranteed 
Time Observation program on 27 March 2004 using all 
four instrument modules. It was re-observed on 23 April 
2005 with the short-low (SL) and long-low (LL) mod- 
ules with increased integration time to achieve a higher 
signal-to-noise ratio (SNR) over the 5 — 36 (im mid-IR 
continuum. A third even longer observation was made on 
16 December 2005 with all four IRS modules. The tar- 
get was acquired using the 22/im (red) peak-up camera in 
high-accuracy mode, and the details of the observations 
are presented in Table [TJ 

The data were processed at the Spitzer Science Cen- 
ter (SSC) (pipeline version 14.0). The two-dimensional 
image data were converted to slope images after lin- 
earization correction, subtraction of darks, and cosmic- 
ray removal. Finally, the data were co-added. In or- 
der to increase the SNR of the subtracted background 
for SL and LL, we combined the background observed 
in off-order and off- nod positions. A detailed explana- 
tion for this method on f aint s ource extraction can be 
found in fWccdman et all ([2006). After subtracting the 
background, the one-dimensional spectra were extracted 
from images with a script version of the Spectral Mod- 
eling, Analysi s and Reduction Tool (SMART, ver. 6.0.4, 
iHigdon et alj (|2004[ )'). We used tapered column extrac- 
tion starting from intermediate pipeline product droop 
files, which only lack stray light and flat-field correction. 
The data from short-high (SH) and long- high (LH) droop 
files used the full slit extraction method from the mean 
of the combined images. We calibrated the flux densities 
by multiplying the extracted spectrum with the relative 
spectral response function (RSRF), which was created 
from the IRS standard stars, HR6348 for SL, HD173511 
for LL and £ Dra for SH and LH, for w hich accurate 
templates are available (|Cohen et al.ll2003l ). 

2.2. Spitzer Imaging with IRAC and MIPS 

The galaxy was imaged at 3.6, 4.5, 5.8 and 8 /im using 
IRAC on 3 April 2004, as well as at 24, 70 and 160 mm 
using MIPS on 7 April 2004 (jEngelbracht et all 12005]) 
(PID:59) (See Table [TJ. The IRAC high dynamic range 
mode was used with a 4 point small cycling pattern of 30 
seconds exposure time for each frame. This resulted in an 
on-source time of 120 seconds for each IRAC filter. The 
MIPS photometry mode was used in small fields with 1 
cycle x 3 seconds at 24 /zm and 2 cycles x 10 seconds at 
both 70 and 160 /im. Two offset positions (±12") were 

7 The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 



TABLE 1 
Spizer Observations of I Zw 18 



AORKEY 


Date 


Instrument 


On-source Time 
(sec) 


9008640 


2004-03-27 


IRS 


84 






IRS (LL) 


120 






TRS fSH~l 


240 






IRS (LH) 


120 


4330759 


2004-04-03 


IRAC (3.6 /im) 


120 






IRAC (4.5 /im) 


120 






IRAC (5.8 /im) 


120 






IRAC (8.0 /im) 


120 


4349184 


2004-04-07 


MIPS (24 /im) 


48 






MIPS (70 /an) 


231 






MIPS (160 /tm) 


42 


12622848 


2005-04-23 


IRS (SL) 


480 






IRS (LL) 


240 


16205568 


2005-12-16 


IRS (SL) 


2040 






IRS (LL) 


840 






IRS (SH) 


2880 






IRS (LH) 


1440 



used to allow proper subtraction of bad pixels. The to- 
tal on-source times were 48, 231 and 42 seconds for the 
MIPS 24, 70 and 160 /im bands. The imaging data were 
processed by the SSC pipeline version 14.0 and the final 
mosaics were obtained from the Spitzer archive. 



3. RESULTS 
3.1. Mid-IR Morphology 

Many ground-based and space-born instruments have 
been use d to obtain optical / UV t o near-IR data for 
IZwl 8 ([Hunter fc Thronsonl 119951: lAloisi et all 119991 ; 
OltTinl l2Tj0ril: iCannon et alJ l200l iHunt et all 120031 
Izotov fe Thuanll2004[ ). iHunter fe Thronsonl (|1995l ) first 
resolved the "main body" of I Zw 18 into stars using the 
Hubble Space Telescope (HST). They have also detected 
filaments of ionized gas up to 450 pc from the center of 
the galaxy. Keck II spectra revealed Ha emis sion as far as 
~1800pc from the "main body" of I Zw 18. Ilzotov et al.l 
(|2001h have also shown that the equivalent widths of 
emission lines are large in this extended envelope. This, 
together with the optical and near-IR colors, suggests 
that ionized gas dominates the emission in the outermost 
regions. 

In Figure 1 we present images of the "mainbody" of 
IZw 18 in four infrared bands (3.6, 4.5, 8 and 24 /im). At 
3.6 /im, where most of the light is due to the stellar pho- 
tosphcric emission, the morphology of I Z w 18 is very sim- 
ilar to that in deep near-IR imaging (IHunt et al.ll2003[) 
and b road band optical optical imaging ([Izotov fe Thuara 
2004). The NW component is noticeably more extended 
and brighter than the SE one. At 8/im the components 
are still clearly resolved while the contrast in the bright- 
ness between the two components has decreased. In nor- 
mal star forming galaxies the emission sampled by the 
IRAC 8 /im filter is typically dominated by dust contin- 
uum, and PAH emission when PAHs are present. Some 
continuum emission from the nearly Rayleigh- Jeans tail 
of stellar photospheric emission may also be present, even 
though its contributio n is typically smal l in late type or 
irregular galaxies (see lSmith et al.ll2007l ). As we discuss 
in the following section, no PAH features are detected in 
the IRS spectrum of IZw 18 down to our 1 a ^0.2 mJy 
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Fig. 1. — IRAC 3.6, 5.8, 8 /im and MIPS 24 /im images of IZwl8 as observed by Spitzer. The galaxy is resolved into a NW and a SE 
component at 3.6 fim which gradually become blended to a single extended source at 24 (im. At the distance of 12.6 Mpc, 1" corresponds 
to ~61 pc. The image is centered on the optical centroid of the galaxy (a = 09°34'02.0", <5 = 55°14'28"). Note the change in morphology 
and peak of the emission as the wavelength increases. 

sensitivity limit. To estimate the contribution of the stel- 
lar continuum to the observed 8 umflux density we follow 
the approach of Uackson et al.l (|2006[ ) and apply a scale 
factor of 0.4 in the 4.5 /zm emission from the galaxy. This 
suggests that no more th an ~25% of the 8 urn fl ux can 
be attributed to the stars (|Engelbracht et aLl l2005) . This 
was also to be expected given the observed steeply ris- 
ing slope of the mid-IR spectrum. Therefore most of the 
"mainbody" emission seen in the 8/xm band is due to 
dust continuum emission. Hence we interpret this grad- 
ual shift in brightness from the NW to the SE compo- 
nent to the probable presence of more embedded star 
formation in the SE component which was obscured in 
the optical broad band imaging. Interestingly, while op- 
tical recombination line ratios give an average extinc- 
tion of only Ay~ 0.2 mag, there are also statistically sig- 



nificantly high Ha/H/3 flux ratios (~3.4, corresponding 
to_Av~ 0.5 mag) in the SE component (|Cannon et al.l 
2002) , suggesting the existence of an appreciable amount 
of dust within the galaxy. Moreover, as can be seen in 
Figure 2, there is some extended 8/im emission, though 
at low levels, to the west of the NW component. This 
emission has similar morphology as the radio contin- 
uum emission detec ted i n the X-band and L-band by 
ICannon et all (|2005f ) and lHunt et~aTl (|2005bl ). which has 
been attributed to low-frequency flux from a synchrotron 
halo. 

In Figure 2, we show contour overlays of the 8 /im and 
24 /im emission on the HST V-band image of the "main 
body" of I Zw 18. At 8 /jm, the source is clearly resolved 
into two components. The centroids of the two compo- 
nents, especially the NW one, are slightly displaced from 
their optical counterparts. Moreover, this displacement 
becomes even more pronounced at longer wavelengths. 
At 24 /im, the two components arc blended into a single 
source, the centroid of which is located slightly closer to 
the SE region. This displacement is real and indicates the 
presence of more 24 /im dust emission in the SE cluster. 
To confirm this change in morphology we convolved the 
8 /im image to the ~ 5.4" size of the 24 /im point-spread- 
function (PSF). Even though the resulting marginally 
resolved source is also elongated in the SE to NW direc- 
tion, the peak emission was further to the NW than the 
peak of the 24/im image. This suggests an actual change 
in the spatial distribution of the various dust temper- 
ature components in the galaxy. In Figure 2, we also 
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Fig. 2. — a) Top: Contour overlay of the IRAC 8 fim image 
on the HST/WFPC2 F555W image taken from the HST archive. 
The astrometric calibration of the HST image has been derived 
using stars from the U.S. Naval Observ atory Astrometric C atalog 
B1.0, and is the same image as used in Hunt et al. (2005bj). The 
contours range from 4<r (<t=0.03 MJysr -1 ) above the sky level 
(1.40 MJysr -1 ) to the peak value of 1.71 MJysr^ 1 . The location of 
the IRS SL (3.6" in width) and LL (10.5") slits, are also indicated 
with the dashed and dash-dotted lines respectively. Note that due 
to the extent of the galaxy, and the fixed position angle of the slit, 
part of both the NW and SE components are not fully covered by 
the SL slit. The size of the PSF at 8 fim (1.8") is also shown at 
the lower right of the panel, b) Bottom: the same optical image 
with the contours of MIPS 24 fira emission (from 4rr and above). 
The size of the PSF is 5.4". 



4 



Wu et al. 




0.000 



15 20 25 30 
Rest Wavelength (/im) 

Fig. 3. — Spitzer/IRS 5-36 /im low-resolution spectrum of 
IZwl8. No scaling factors have been applied to stitch the dif- 
ferent orders and modules. We indicate several of the well known 
mid-IR fine structure emission lines, detected in the spectrum (see 
also Fig. 4). 

overlay the IRS SL/LL slit on the image of the galaxy. 
As one can see, since the SL slit is only 3.6" wide, some of 
the flux from both components is not properly sampled, 
resulting in an underestimate of the extended emission 
from the galaxy. However, because the spectrograph is 
not sensitive to the low-surface brightness emission from 
the areas denoted by the lowest 8/xm contours, this does 
not bias our analysis of the global spectral properties of 
the system or affects any of the conclusions drawn. 

3.2. Mid-IR Spectral Features 

Figure 3 shows the 5.3-36 /im low-resolution spectrum 
of I Zw 18 as observed by t he IRS. The SNR is ~3 times 
higher than that shown bv lWu et alj (|2006f) . As we will 
discuss in §3.4, the global shape of the mid-IR spectrum, 
reveals that the IRS spectrum of I Zw 18 continues to rise 
steeply with wavelength from 5 /im all the way to 36 /im. 
This i s unlike the case of SBS 0335-052E ([Houck et all 
2004b), the third lowest metallicity galaxy to date, which 
has a nearly flat continuum peaking at ~28 /im in f„. 

The improved SNR of the new spectrum enables us 
to detect for the first time several mid-IR forbidden 
lines. Fine structure lines, such as [SIV] at 10.51 /im and 
[Nelll] at 15.55 /im can clearly be seen, even in the low- 
resolution spectrum. Several additional forbidden lines, 
such as [Nell] at 12.81 /im, [SIII] at 18.71 and 33.48 /im, 
as well as [Sill] at 34.82 /im are identified in the high- 
resolution spectrum (see Figure 4). [OIV] at 25.89 /im 
and [Fell] at 25.99 /im are blended in LL, but can clearly 
be resolved in LH. The line fluxes measured from the IRS 
high-resolution spectrum are reported in Table [5] The 
observed line ratio of [NcIII]/[NeII] is ~ 5 and the ratio 
of [SIV] /[SIII] (18.71 /im) is - 2. This indicates that the 
radiation field in I Zw 18 is much harder than in a typi- 
cal starburst galaxy, where [NeIII]/[NcII] is usually < 1 



galaxy, 

(|Brandl et alj [20061 ) . and eve n harder than the majority 
of the BCDs observed so far (|Hunt et al.ll2006HWu efaLl 
[20061) . 

As was the case with SBS 0335-052E, the 5-15 /im spec- 
trum of IZwl8 does not show any detectable Poly- 
cyclic Aromatic Hydrocarbon (PAH) emission. Using the 
new mid-IR spectrum, we measure a 3<r upper limit of 
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Fig. 4.— Mid-IR fine structure lines of [SIV] (10.51/im), 
[Nell] (12.81/im), [Nelll] (15.55/im), [SIII] (18.71, 33.42/^m) [OIV] 
(25.89/im), [Fell] (25.99/im) and [Sill] (34.82/tm) from the high- 
resolution spectrum of I Zw 18. Note that the sky emission has not 
been subtracted. 

3.5xlO~ 22 Wcm~ 2 for the 6.2 /mi PAH feature, as well 
as an equivalent width (EW) of <0.23/im. The 11.2 /im 
PAH emission has an upper limit of 1.3x 10 _22 Wcm -2 
and EW of <0.12/im. The decreasing strength of 
the PAH features in low met allicity environment is a 
current topic o f interest (see [Engelbracht et all 120051: 
Wu et alj I2006t lO'Halloran et alj \200Q. Uackson et al l 
20061: iMadden et alj|2006t iBeirao et al.ll2006l) . The exact 
reason why PAHs tend to be absent in low-metallicity, 
high-excitation star-forming regions is not yet clear, al- 
though it almost certainly has to do with a combination 
of several effects, such as low carbon abundance, shock 
destruction of grains by supernovae, high ionization or 
excitation resulting from low metallicity and/or extreme 
radiation field intensity. 

3.3. Neon and Sulfur Abundances 

The fine-structure lines of sulfur and neon measured 
in the high-resolution spectrum can be used to derive 
the ionic abundances relative to hydrogen. To derive 
these abundances one needs at least one hydrogen re- 
combination line (usually H/3), as well as an estimate 
of the electron temperature (T e ) and electron density 
(N e ). The H/3 fluxes for IZwl8 reported in the liter- 
ature often correspond to small regions of the galaxy 
or of one of the two components. The IRS SH slit 
from which the IR lines were measured contains part of 
both components. We derive the H/3 flux from the ther- 
mal c omponent of the 3.6 cm continuum (|Cannon et al.l 
2005), which is not affected by extinction effects and 
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encompasses the entire "main body". The SH slit in- 
cludes ~52% of the total H/3 emission and this gives an 
H/3 flux of 6.1xl0 _14 ergscm~ 2 s~ 1 . The reported elec- 
tron temperatures in the lit erature range from 18 000 to 
20000 K for O II and HI (ISkillman fc Kennicuttlll993l: 
Izotqv fc Thuanl 119991 : iThuan fc Izotovi 120051 : IShi et all 
2005J), with "slightly lower value for S II (T e ~15 500 
K). The infrared lines are less sensitive to T e and here 
we adopt an average temperature of 19 000 K. An elec- 
tron density (N e ) of 100 cm " 3 was assumed (|Shi et al.l 
l2005t llzotov fc Thuan| [l999?). Using the above quantities 
and resolving the population of the l evels for each ion, the 
ionic a bundances can be derived (see lBernard Salas et al.1 
(120011 ). their eq. 1). 

The calculated ionic abundances are presented in Ta- 
ble [21 The elemental abundances can be obtained by 
adding the contribution of each ion. In the case of neon, 
only Ne II and Ne III are observed. The presence of 
O IV suggests that some Ne IV may be present. As- 
suming that the contribution of Ne IV is 10% of that 
of Ne HI as has been found in s ome planetary nebulae 
(PN 8 . lBernard-Salas et all (120031 )). the total neon abun- 
dance is 5.3xl0 -6 , which is 1/23 Z . Similarly, adding 
up the ionic abundance of S III and S IV, the total sul- 
fur abundance is 4.0xl0 -7 , which is 1/35 Z Q . How- 
ever, this might be a lower limit because sulfur abun- 
dance as derived from galactic planetary nebulae and 
Hli regions is lower than the s olar sulfur abundanc e 
(iPottasch fc Bernard-SaTasll200fl iMaciel fc Costal 120031 ). 
probably due to an overestimate of the solar value. Com- 
paring directly with the ionic abundance from the optical 
(llzotov fc Thuanl [19991. our Ne III abundance is nearly 
twice as high while S III abundance agrees quite well. 
The difference in Ne III could come from the electron 
temperature that is used to derived the ionic abundance. 
The optical is more affected by the change in temperature 
as compared to the infrared. Lowering T e from 19000 K 
to 15000 K would double the Ne III abundance derived 
from the 3869A line. It is known that in some PN, the 
temperature obtain ed from the Ne HI ion is l ower than 
that from the O III ([Bernard Salas et al.l l2002). Another 
possible explanation is that there are some regions with 
dust obscuring the optical emission lines. Overall, the 
neon (1/23 Z Q ) and sulfur (1/35 Z Q ) abundances we de- 
rive using the infrared lines are consistent with the nebu- 
lar oxygen abundance (1/30 Z Q ), which supports the low 
extinction (Ay— 0.2 mag) derived from hydrogen recom- 
bination lines bv I Cannon eTai] (I2002T ). 

3.4. Comparison with SBS 0335-052E and NGC7714 

I Zw 18 and SBS 0335-052E share some properties but 
are very different in other aspects. Perhaps the most 
salient difference between their spectral energy distri- 
bution (SED) is the fraction of their luminosities emit- 
ted in the IR. Using the 15 /im and 30 /im flux den- 
sities of I Zw 18 and applying an empir ical relation 
in starburst galaxies (jBrandl et al.1 120061 ). we derive 
L m = 1.8xlO 7 M . While SBS0335-052E has L m ~ 
10 9 Lq and L m /L B - 1.3, IZwl8 has L m - 10 7 
Lq and Lir/Ls ~ 0.3; the relative infrared luminos- 
ity is a factor of 4 times greater in SBS 0335-052E. 

8 The high excitation and densities in BCDs are a closer match to 
PN than the typical H II regions found in normal starburst galaxies. 



TABLE 2 
Fluxes and ionic abundances 3 - 



Arest 


Feature 


Obs. Flux b 


N iott /N p c 


10.51 


[S IV] 


4.8±0.3 


1.3xl0~ 7 


12.81 


[Ne II] 


0.9±0.1 


1.4xl0" 6 


15.55 


[Ne III] 


4.6±0.2 


3.6xl0" 6 


18.71 d 


[S III] 


2.3±0.2 


2.8xl0" 7 


25.89 


[O IV] 


4.9±0.3 




25.99 


[Fe II] 


3.4±0.3 




33.48 d 


[s mi 


12. Oil. 2 




34.82 


[Si nf 


15.8±1.8 





a Only the abundances of ions used to calculate the metal- 
licities in this work are listed in this tabled In units of 
10 -15 erg cm~ 2 s _1 . c Ionic abundance relative to hydrogen. 1 ^ Here 
we use the [S III] at 18.71 fim when calculating the sulfur abun- 
dance because this line is in the SH module and provides us a more 
reliable measurement as compared to the 33.48 /im line which is 
at the edge of the LH module and is noisy. 

Another important difference in the infrared SEDs of 
IZwl8 and SBS0335-052E is the peak wavelength of 
the SED. Th e SED of SBS 0335- 052E peaks at -28 /im 
in f v space (jHouck et al.l l2004bT) . indicating very little 
cold dust. IZwl8 has a clear detection of 34±2.4mJy 
at 70 /im(Engelbracht et al. 2007, in preparation), and 
the ratio of f7o/f24 is more than a factor of 5 while the 
same ratio in SBS 0335-052 is less than 1. This suggests 
that contrary to SBS 0335-052E, which is a similarly low 
metallicity BCD, there is a significant amount of cold 
dust in I Zw 18. 

We can further compare the properties of I Zw 18 with 
higher-luminosity and more metal-rich starbursts. In 
Figure 5, we present the spectra of I Zw 18 and SBS 0335- 
052E, together with a typical starburst galaxy, NGC7714 
(|Brandl et al .1 [2004h . all of which have been normalized 
to the flux density of I Zw 18 at 22 /im. Setting aside 
the strong PAH emission features present in NGC7714, 
there is a striking similarity between the mid-IR contin- 
uum of IZwl8 and NGC7714, although the latter has 
a metallicity of more than half solar. At short mid-IR 
wavelengths (A <10/im), the warm dust component dom- 
inates and it does not change much ev en for an extended 
starburst galaxy ([Brandl et al.l l2006h . At longer wave- 
lengths, emission from larger cooler grains is dominant. 
After normalization at 22 /im, the 70 /im flux of I Zw 18 
differs less than 20% when compared to the 70 /im flux 
of 9.53 Jy (Engelbracht, private communication) or the 
IRAS 60 /im flux of 11.16 Jy for NGC7714. Figure 5 
suggests that similarly low metallicity galaxies can have 
both very flat or very steep spectral slope, while in the 
latter case, the spectral slope can be as steep as that of a 
typical starburst. This leads us to conclude that metal- 
licity is not the main parameter driving the difference in 
the shape of the mid-IR spectral slope (|Wu et al.ll2006h 
and the infrared part of the SEDs. 

Based on the new Spitzer measurements of the mid- 
and far-infrared emission from I Zw 18 one could in prin- 
ciple attempt to model the global SED of the galaxy. 
We did e xplore this avenue using modelling tools such as 
DUSTY (jlvezic et al.lll999l ). but with only limited suc- 
cess. The complex geometry of IZwl8 and the large num- 
ber of free parameters in the available models prevented 
us from significantly constraining the physical conditions 
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Fig. 5.— The mid-IR spectra of IZwl8, SBS 0335-052E and 
NGC7714 between 5-36 fim. The spectra have been normalized to 
the flux density of I Zw 18 at 22 /an. SBS 0335-052E stands out 
as having a very fiat continuum in the mid-IR, indicating there is 
much less cold dust in this galaxy. I Zw 18 and NGC7714 display 
a striking resemblance in this wavelength range. 

of the dust. We thus refrain from elaborating on these 
results until more data are available. 

4. STAR FORMATION RATE IN I ZW 18 

Deriving the star formation rate in nearby galaxies 
from various observational indicators and understanding 
the possible variations in the results is extremely useful 
for "predicting" the properties of high-redshift galaxy 
populations, where only sparse data are available. Us- 
ing the available data, we calculated the SFRs of I Zw 18 
from different indicators and we present our results in 
Table El 

The SFR estimated from t he Ha luminosit y gives 
a value of 0.05 M^yr" 1 (jKennicutt et all 119941 : 
ICannon et al.l 120021 ). However, IZwl8 has only 
1/30 Z Q and lower metallicities may result in a re- 
duced SFR for a given Ha luminosity (|Lee et al.l l2002t 
iRosenberg et alj l2006f). Using the metallicity correc- 
tion recipe of iLee et al.1 (|2002| ). the SFR would be 
-0.03 Moyr" 1 . 

The radio continuum emission is another important di- 
agnostic of star-formation process and it is not affected 
by dust extinction effects. The thermal free-free emis- 
sion is a direct indicator of SFR and it is typically only 
— 10% of the total radio con tinuum at 1.4 GHz for nor- 
mal galaxies (jCondonl fl992). However, in IZwl8, the 
fraction of ther mal component is three times the typica l 
value (-30%: iHunt et all l2005bl : ICannon et all l200l . 
Using the relation between radio fre e-free emission an d 
SFR, we derive a SFR=0.1 Moyr" 1 fHunt et al.ll2005bh . 
The non-thermal component of the radio continuum can 
also be used t o calculate SF R and we derive a SFR of 
0.03 Mgyr" 1 (jCondonl IT992I ). a factor of 3 lower than 
the thermally-derived SFR. This difference is a direct 
consequence of the unusual value of the "thermal/non- 
thermal" fraction in I Zw 18. This is an important caveat 
that should be considered when applying the standard 
correlations, which have been established for normal star 
forming galaxies, in very young low- metallicity systems. 

Finally, the SFRs estimated from the infrared are 
significantly lower. Using the total infrared luminos- 
ity of 1.8xl0 7 L Q , we derive a SFR=0.003M Q yr- 1 



TABLE 3 
Star Formation Rate Estimates of 
IZwl8 



SFR indicator SFR Reference 
(Moyr" 1 ) 



L Ha b 


0.05 


(1) 


L Ha c 


0.03 


(2) 


Ltir 


0.003 


(3) 


L24^m 


0.006 


(1) 


Lir+uv 


0.02 


(5) 


Bt hernial 


0.1 


(6) 


Lnon — thermal 


0.03 


(6) 



a References : 1) IKennicut t_et_alj l|1994f) , 
2) ILee et alj <|2002D 3 ) IKennicut3 1199811 
4) IWu et alj 1 120051 ) 5) IBell et alj J2005D 6) 
ICondonl U992l) . b Not corrected for low- 
mctallicity effects. c After correcting the low- 
mctallicity effects. 

(|Kennicuttl I1998T) . while using t he 24 /zm em i ssion , 
we find that S FR=0.006M w yr- 1 (|Calzetti et al.1 12001 
IWu et all |2005| ). This is probably because the dust 
content in I Zw 18 is so low while the above relations 
have been calibrated for sources of high optical depth, 
where virtually all of the UV radiation is converted to 
infrared luminosity. A simple calculation using the red- 
dening curve of the Small Magellanic Cloud assuming an 
Ay =0.2 mag suggests that a significant amount of UV 
light has leaked out without being absorbed by the dust, 
thus the lower SFRs estimated from the infrared are not 
unexpected. The readers should be aware of these com- 
plications when applying the canonical infrared relations 
for estimating the SFR in environments with low dust 
optical depth. If we were to assume that Lj R accounts 
for the bolomet ric luminosity of th e obscured popula- 
tions while Luv (jKinnev et al.lll993t ) represents the con- 
tribution o f the unobscured stars, and use the eq.l of 
IBell et al.1 (|2005h 9 . we find a SFR of 0.02 Moy^ 1 , more 
consistent with the SFRs derived using the Ha or radio 
luminosities. 

5. CONCLUSIONS 

We have explored the mid-IR and far-IR properties of 
the archetype BCD I Zw 18 based on Spitzer data: 

1. Using the low-resolution modules of the IRS, we have 
acquired the deepest mid-IR spectrum of this galaxy ob- 
tained so far. No PAH emission is found which confirms 
the absence of PAHs in low metallicity systems. How- 
ever, the mid- to far-IR spectral slope of I Zw 18 is sur- 
prisingly similar to NGC7714, a typical starburst galaxy 
with half solar metallicity. This, especially the MIPS 
70 /im detection, would suggest the presence of a signifi- 
cant amount of cold dust in IZwl8. 

2. Variations in the morphology of the galaxy from 
3.6 to 24/zm imaging imply that more dust emission is 
present in its SE component t han in the NW one. This 
agrees well with the results of ICannon et al.l (|2002h . 

3. The mid-IR fine-structure lines identified in the 
high-resolution spectrum of I Zw 18 imply a neon and 
sulfur abundance of 1/23 and 1/35 Zq respectively, con- 
sistent with the optically derived oxygen abundance of 
1/30Z Q . 

9 SFR=9.8xlO~ 11 (L IR +2.2 L uv ) 



Spitzer view of I Zw 18 
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4. Estimates of the star formation rates calculated 
from different indicators show considerable scatter. Lir 
and L24 Aim give lower SFRs when compared with results 
using Ha or Li.4qhz, probably because the low dust con- 
tent in this galaxy can only convert a small fraction of 
the UV radiation emitted by stars into Lir. This should 
be considered when interpreting star formation rates de- 
rived for high-redshift low metallicity galaxies. 



We thank Chad Engelbracht for graciously provid- 
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our source before publishing the data. We also thank 
Rob Kennicutt for insightful discussions. We thank the 
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for this work was provided by NASA through Contract 
Number 1257184 issued by JPL/Caltech. 
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